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Abstract 

We propose the first algebraic determinantal formula to enumerate tilings of a 
centro-symmetric octagon of any size by rhombi. This result uses the Gessel-Viennot 
technique and generalizes to any octagon a formula given by Elnitsky in a special 
case. 
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Introduction 



The enumeration of tilings of a centro-symmetric polygon by rhombi is a no- 
toriously difficult problem that concerns discrete mathematics and theoretical 
computer science, as well as theoretical physics, in relation with quasicrystal- 
lography. In the latter community, these tilings are usually called "random 
tilings with octagonal symmetry". We address the following issue: given a 
centro-symmetric octagon Oa,b,c,d, of integral sides lengths a, b, c and d (read 
clockwise; see figure 2, left), in how many ways is it possible to fill it entirely, 
without any gap or overlap, with the following six species of tiles: two dif- 
ferently oriented squares, and four differently oriented 45° rhombi, the six 
of them with unitary side lengths? So far, this question has been solved in 
very particular instances only. We denote by Ta^b,c,d the set of all the tilings of 
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Oa,b,c,d: and by TaM,c,d = \Xi,b,c,d\ the cardinality of Ta,b,c,d- For example, figure 1 
displays the eight tilings of the set 







Fig. 1. The eight tihngs of the set T^, 1,1,1. 

Small systems have been studied in references [1,2] up to sizes of some hundred 
tiles (see table 1). However, the technique employed cannot reasonably provide 
tiling enumerations for bigger octagons. On the other hand, Elnitsky gave in 

Table 1 

Some tihng enumerations computed in ref. [2]. The number of rhombi is given in 
column 3. 



a, b, c, d 


'^a^byC^d 


# tiles 


1,1,1,1 


8 (see fig. 1) 


6 




5383 


24 


3,3,3,3 


273976272 


54 


4,4,4,4 


1043065776718923 


96 


5,5,5,5 


296755610108278480324496 


150 



ref. [3] two formulas when two sides of the octagon are set to 1: 



and 



T 

-^a,b,l, 



2(a + 6 + 1)! (a + 6 + 2)! 
a! 6! (a + 2)! (6 + 2)! ' 



(2) 



The first formula has been later partially simplified [2] : 



7^ 

-'- a,l,c,l 



(a + c+1)! 



a! c! (2a+l)(2c+l) 



2(a + c+l)! 



ai c! 



+ E 

k=0 



a\ c 



2k - 1 \k \k, 



(3) 
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where the last sum can be written in terms of a hypergeometric function 
^ 3F2[-l/2,-a,-c;l/2,l;l]. 




fe=0 



2k-l \k \k, 



(4) 



We propose a generahzation of the first formula (1) to any side lengths, where 
Ta,b,c,d is written as a sum of products of determinants. Even if the complexity 
of our formula increases with the system size, it is the first explicit algebraic 
expression to count tilings of an octagon (see eq. (12)), which can be in prin- 
ciple calculated for any system size. 

As it is discussed below into detail, Elnitsky's proof uses a "square grid rep- 
resentation" of tilings (see figure 2), which is closely related to the "de Bruijn 
dualization" , a wide-spread technique in quasicrystal science. This dualization 
has proved powerful to handle rhombus tilings in several circumstances. The 
present paper also uses this technique, thus generalizing Elnitsky's proof. 




Fig. 2. Left: Example of octagonal tiling of a centro-symmetric octagon of sides 
a, b, c, d. There are 6 species of tiles with unitary side length: two squares and four 
45° rhombi. One species of squares, the "tilted" squares, is emphasized in medium 
gray. They lie at the intersections of the light gray and dark gray de Bruijn lines. De 
Bruijn lines are defined in section 1 and are made up of adjacent rhombi sharing an 
edge with a given orientation. Right: square grid representation of the same tiling, 
obtained by shrinking all colored tiles; see section 1. The white disks on the grid 
keep track of the position of the tilted squares of the original tiling. The orientation 
of edges (arrows) will be discussed later. 

We now state our main result. Given the side lengths a, b, c and d, we denote 
by X (resp. Y) the set of families of integers {xk,i) (resp. {yk,i)), k = 1, . . . ,b, 
I = 1, . . . ,d, satisfying the relations: 



< Xk,i < a, 

Xk,i ^ Xk'^v if k < k' and / < 
< yk,i < c, 

yk,i < yk',i' if k > k' and / < 



(5) 
(6) 
(7) 



Note that conditions (6) and (8) are not exactly similar. In the following, these 



3 



integers will be the coordinates of the white disks in the grid representation 
(see figure 2, right, and section 1). In addition, we set by convention for k — 
1, . . . ,b and I = 1, . . . ,d 

Xkfl = 0, Vkfl ^ 0, 
Xo,i = 0, ?/o,/ = c, 
Xb+i^i = a, 7/6+1,/ = 0, 

Xk,d+1 — 0.1 yk,d+l — C- 



The reasons for this convention will be explicited below. For any two such 
sequences x = {xk,i) and y = {yk,i), we define the matrices M^'^^(x,y) and 
P^'"^ {x, y) as follows: M*^") (x, is a 6 x 6 matrix of coefficients 

m,^ = h^"" ~ ^^'"-1 + 7 (10) 



for l<i,j <h and P^''\x,y) is a, d X d matrix of coefficients 

— -I- ■ — ■ J ' \^^) 



for 1 < i,j < d. Note that, by convention, we set (^^^ — whenever ^4 < or 
B <0or B > A. 

Theorem 1 With the previous notations, the number of rhombus tilings of a 
centro- symmetric octagon of sides a,b,c,d reads: 

Ta,b,c,d= E ndetM(")(a;,|/)ndetP('-'Ha;,i/). (12) 

{x,y)eXxY u=l v=l 



It is demonstrated below that the determinants come from the enumeration, 
by the Gessel-Viennot method (presented below), of tilings of independent 
sub-domains of the octagon delimited by some points of coordinates {xk,i, yk,i) 
in the square grid representation. 

When b = d — 1, the previous expression is reduced to Elsnitsky's relation (1). 

Note that by contrast, relation (2) is not a spacial case of this formula. Beyond 
this simple case, the number of terms in the formula grows with the octagon 
size. For example, for (a, b, c, d) = (2, 2, 2, 1), the formula contains 6 x 6 = 36 
terms to count the 480 tilings. For {a,b,c,d) = (2,2,2,2), there are 20 x 
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20 = 400 terms and 5383 tilings. More generally, the number of terms grows 
exponentially with the number of tiles, but it nevertheless grows exponentially 
more slowly than the number of tilings. As a consequence, this formula is 
exponentially more compact than the crude enumeration of tilings. This point 
is discussed in the conclusion. 



1 Octagonal tilings and the square grid representation 

In this section, we show that octagonal tilings are conveniently represented by 
families of directed paths running on a rectangular patch of square grid. This 
representation was used by Elnitsky [3] and it is reminiscent of the prior "de 
Bruijn dualization" [4,5,6,7] and derived representations [1,2]. We first expose 
briefly the de Bruijn dualization process. Figures 2 and 3 will help the reader. 
To begin with, we notice that tile edges can have four possible orientations. 
We define a family of de Bruijn lines for each orientation: de Bruijn lines are 
made up of adjacent rhombi sharing an edge with this orientation. It is always 
possible to extend these lines through the whole tiling up to a boundary tile. 
Two examples of lines are presented in both figures 2 and 3, belonging to two 
different families. A rhombic tile is situated at the intersection of two lines 
of different families (see the figures). By construction, lines of a same family 
never intersect; there are respectively a, b, c and d lines in each family. 

Now we show how to translate the de Bruijn's representation of a tiling into 
its square grid representation. In figure 3, we show this correspondence in the 
simplest case b = d = 1 [3]. The idea is to shrink the de Bruijn lines of two 
families among four, so that they become paths on a square grid, as displayed 
in figures 2 or 3. Because all tiles of a de Bruijn line have an edge of a given 
orientation, these paths are directed. The b paths of the first family (denoted 
by SW) go from the south-west corner to the north-east one (dark gray); 
they can follow eastward and northward edges only; the d paths of the second 
family (denoted by NW) go from the north-west corner to the south-east one 
(light gray); they can follow eastward and southward edges only. 

In the simplest case b = d = 1, to avoid ambiguity due to path tangency and to 
make this correspondence bijective [3], we keep track of the intersection of the 
de Bruijn lines thanks to a distinguished vertex, represented by a white disk 
in the right figure. It marks the position of the unique tilted square (medium 
gray). When 6>lord>lasin figure 2, there are bd intersections and 
therefore bd tilted squares. Each of them must be located by a distinguished 
vertex on the square grid. Paths do not cross in a same family even though they 
can be locally adjacent (see figure 2). We denote the paths of SW (resp. NW) 
by SWi, . . . , SWb (resp. NWi, . . . , NWd) from left to right. As a consequence, 
distinguished vertices are indexed by two integers k and I, and are denoted by 
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Fig. 3. Example of tiling in the case h = d = \ (left) together with its square 
grid representation (right). We have emphasized two de Bruijn lines of the tiling, 
belonging to two different families. Their intersection coincides with the gray tilted 
square. They are represented by two directed paths joining two opposite corners 
of the grid. These paths can have multiple intersections. To avoid ambiguity on 
the position of the gray tilted square on the original tiling, a white circle marks a 
distinguished vertex to keep track of its position on the grid. Note that when two 
paths are tangent, we have slightly shifted one of them for the sake of readability. 
However, these paths are in fact superimposed and lie on the same grid edge. 

DVk,i (see figure 4). 



2 The Gessel-Viennot method 



The Gessel-Viennot method [8,9] is a combinatorial technique for the count- 
ing of configurations of directed non-intersecting paths on oriented graphs. 
This technique has already proved very useful for the enumeration of rhombus 
tilings (see [3,10] for examples, as well as section 3 of the present paper). It 
has been extensively described in the literature [8,9] and we shall only briefly 
explain it in the present paper, focusing on the underlying ideas and not on 
technical details. The method is rather general and can be applied to any 
acyclic oriented graph in which are selected two families of vertices, di 
("departure" vertices) and aj ("arrival" vertices), i, j = 1, . . . ,n. We consider 
directed paths, running on Q, starting from one vertex di and arriving at one 
vertex aj. By "directed", we naturally mean that the paths must follow the 
edge orientations. In addition, this graph is supposed to satisfy the property 
of compatibility, if two directed paths on Q are going respectively from to 
ttji and from di^ to Uj^, if these paths do not cross, and if ii < 12 then ji < j2- 
This property is very specific to two-dimensional graphs. 

We are interested in the number Dn of configurations of n non-intersecting 
directed paths on Q, where the i-th path goes from di to Oj: two paths are 
said to be non-intersecting if they share no vertex; n paths are said to be 
non-intersecting if any two paths are non-intersecting. If we denote by Xij the 
number of paths going from di to aj, then the Gessel-Viennot theorem states 
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that 

Dn = det(Aij) 



(13) 



The idea of the proof is that in this determinant, aU configurations of n paths, 
whether intersecting or not, the i-th path going from di to a^^^^, for any per- 
mutation (7, are counted, with a + or — sign. Because of these signs, all 
configurations with one or more intersections cancel two by two. Only the 
non-intersecting configurations remain. They are exactly the configurations 
under interest thanks to the property of compatibility. The interested reader 
is referred to Stembridge [9] for more detailed explanations. 



3 Proof of theorem 1 

We are now ready to prove theorem 1. First of all, we need to endow the square 
grid with integer coordinates in order to locate the positions of distinguished 
vertices. They are defined according to the usual conventions, so that the 
south-west and north-east corners have respective coordinates (0, 0) and (a, c). 
The coordinates oiDVk,i are denoted by {xk,h yk,i)- By extension, we also define 
the vertices DVkfi = (0,0), DVk,d+i = (a, c), DVq^i = (0, c) and -DH+i,; = 
(a, 0). They are the ends of paths of families SW and NW . If their coordinates 
are also denoted by {xk,i,yk,i), these last definitions are compatible with the 
conventions (9) given in introduction. All these coordinates naturally obey 
relations (5) and (7). 

Furthermore, because of the directed character of de Bruijn lines, distinguished 
vertices DVkj are constrained by some conditions when they belong to the 
same paths, and they must obey relations (6) and (8) as well. These four 
conditions define the sets X and Y, as it was stated in the introductory part. 

Let X = (xk^i) G X and y = {yk,i) G F be an admissible set of coordinates 
of the distinguished vertices. We denote by %:^y the subset of tilings of Ta^b,c,d 
in the square grid representations of which the distinguished vertices have 
these coordinates. The subsets are two-by-two disjoint so that Tafi,c,d — 
J2{x,y)€XxY \%:,y\- Our purposc is now to calculate each \%:^y\. This calculation 
is feasible because for a given {x,y), the subset %;^y can be factorized into 
simple sets (see eq. (15)). Each of them can in turn be counted by the Gessel- 
Viennot method, which leads to relation (12). 

First we need to introduce two definitions. Given a configuration of vertices 
DVk^i = {xk,i,yk,i), we fix / and we consider in isolation the vertices DVk^i-i as 
well as DVk,u ^ = 1, • • • , & (see figures 4 and 5). Then we define the set sw{l) 
of all the configurations of b directed non-crossing paths, the k-th path going 
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Fig. 4. Examples of configurations of non-crossing directed paths, belonging to the 
sets sw{l) (left) and sw{l) (right). Each path (solid black lines) is a section of a path 
SWk and goes from DVk^i^i to DVk^i (white circles). By convention, DVkfi = (0,0). 
The gray paths NWi, NWi^i and NWi do not belong to the path configurations 
and are drawn for information only. The dashed black lines represent the possible 
continuations of the original paths SW^- In this figure also, we have slightly shifted 
tangent paths for the sake of readability, whereas in reality they run on the same 
grid edge. 

from DVk,i-i to DVk^i, with k = 1, . . . ,b. These paths are directed from south- 
west to north-east. They have no constraint except that they are directed and 
non-crossing (these paths can have tangencies). In a similar way, we define 
the sets nw{k) for any k: they are the sets of all configurations of d directed 
non-crossing paths, going from north-west to south-east. The l-th path goes 
from DVk-i,i to DVk,i. 

Now in order to prove theorem 1, we start from the following observation, 
illustrated in figure 4: in the distinguished vertices DVk^i are held fixed 
and one can consider independent patches of the families SW or NW, as 
follows. Without loss of generality, we focus on SW. We cut each path SWk 
into d + 1 sections, denoted by SWk{l), where / = 1, . . . , ci + 1. The section 
SWk{l) goes from DVk^i-i to DVk^i- Then all the h sections SWkil), k = 1, . . . ,b 
form a local path configuration denoted by Psw{l)- It belongs to sw{l). In a 
similar way, the corresponding grid patches defined with respect to the family 
NW are denoted by Pj^w{k) and belong to nw{k). Therefore, when x and y 
are fixed, we have the natural inclusion 

d+l b+l 

C n sw^u) W nw{v), (14) 

where the products are direct. We prove below that 
Lemma 1 The previous inclusion is an equality: 

d+l b+l 

^:?/ = n sw{u) Jl nw{v). (15) 

U=l 11=1 
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It follows from eq. (15) that 

disjoint d+1 b+1 

Ta,c,b,d= u n ^'^(^) n ^^(^)' (16) 

{x,y)^XxY «=1 v=l 



and that 

d+l b+1 
[x,y)eXxY u=l v=l 



The remainder of the proof consists in calculating the cardinalities [^^(li)! 
and \nw{v)\ by the Gessel-Viennot method. Indeed, it is also demonstrated 
below that 

Lemma 2 When x and y are fixed, 

\sw{u)\ ^ det M^''\x,y) ; \nw {v)\ ^ det P^''\x,y). (18) 

Proof of lemma 1: 

We need to prove the reverse inclusion 

d+l b+1 

%,yZ^'[{sw{u)^nw{v). (19) 

u=l v=l 

Configurations from the sets sw{u) provide sections of paths from DVk^i-i 
to DVk^i- When concatenated, these sections provide complete directed non- 
crossing paths from (0,0) to (a,c), which form a family SW . In a similar 
way, sections from the nw{v) provide directed non-crossing paths from (0, c) 
to (a,0), forming a family NW . We only need to check that any two paths 
from SW and NW only cross at the distinguished vertices DVk^i- This point 
is ensured by the directed character of path sections (see figure 5). This last 
observation is crucial and all the demonstration relies on it: it ensures the 
reverse inclusion and therefore the direct character of the product (15), from 
which our enumerating formula ensues. 



Proof of lemma 2: 



So far we have used the terminologies "non-mtersecitm^' paths" in section 2 and 
"non- cro5sm^ paths" in section 3. Now it is time to precise what aspects these 
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Fig. 5. Example of path concatenation involved in the proof of lemma 1. The first 
line displays configurations from the sets sw{u) (black paths). The second line from 
the sets nw{v) (gray paths). Altogether, the superposition of the six patches forms 
a square grid representation. Black paths and gray paths can cross only at the 
distinguished vertices (white circles). 

two terms cover. We have seen that non-crossing paths can have tangencies, 
that is to say they can share vertices or edges of the grid, but they cannot step 
over one another. In particular, non-crossing paths of families SW (or NW) 
share their ends, but can be indexed from west to east without ambiguity. 

On the contrary, non-intersecting paths cannot share any vertex or edge. 
Therefore, if we want to use the Gessel-Viennot method, we need to transform 
configurations of non-crossing paths on the square grid into configurations of 
non-intersecting paths. The trick consists in shifting non-crossing paths, as 
it is illustrated in figure 6. The trick is standard and was already used by 
Elnitsky [3] for example. 




Fig. 6. Left: a configuration of non- crossing paths; Right: a configuration of 
non-intersecting paths. The second configuration is obtained from the first one via 
the shifts defined in the text. This correspondence is one-to-one. 



We use a unitary shift vector u = (1,-1) and we shift the k-th path section 
SWk{l) by a vector {k — l)u (see the figure). The new paths still belong to 
the square grid. The k-th new path section goes from the new vertex DV^ i_i 
to the new vertex DV^i, of coordinates 

x'k,i = ^k,i + {k - I) and y',^^i = yk,i - {k - 1). (20) 
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The so-obtained path configuration is non-intersecting by construction. This 
correspondence between configurations of non-crossing paths and configura- 
tions of non-intersecting paths is bijective. 

Now we use the Gessel-Viennot technique by setting the departure vertices 
di = DV/i_i and the arrival ones aj = DV-i- The number of directed paths on 
the square grid going from a vertex D = {xd, Vd) to a vertex A = {xa, Ha) is 
simply given by the binomial coefficient 

^ _ (xa- xd + Va-Vd 
\ xa-xb 

Then one computes the coefficients Xij involved in the Gessel-Viennot method: 
\ij = rriij. Thus one obtains the matrices M*^"^ (x, y), the determinants of which 
count the elements of sw{u). In the same way, to count the elements of the sets 
nw{v), one must shift the sections of paths NWi by (/ — l)v where v = (1, 1). 
One gets the matrices P^'"\x, y) and takes their determinant, which completes 
the proof of lemma 2. 



4 Conclusion 

We have demonstrated how Elnitsky's technique can be generalized to oc- 
tagons of any size, leading to an explicit enumerative formula (theorem 1). 

We also notice that the conditions (5) and (6) (resp. (7) and (8)) that define 
the set X (resp. Y) are identical to the conditions defining plane partitions 
of height a (resp. c) on a 6 x d grid [11]. This point is remarkable because 
such plane partitions are known to be equivalent to rhombus tilings filling a 
centro-symmetric hexagon of sides lengths 6, d and c (resp. 6, d and a) [12]. 
We have derived a partial combinatorial interpretation of our formula (12) in 
terms of these tilings of hexagons. It is related to a natural decomposition of 
the configuration sets of tilings of octagons, as described in ref . [2] . But it goes 
beyond the scope of the present paper and will be described elsewhere [13] . 

If T^l^p denotes the number of tihngs of the centro-symmetric hexagon of sides 
a, h and c, the previous remark leads to the lower bound 

rp ^ fT-ihex ^hex 

J-a,b,c,d ^ ^b,d,c ^b,d,a-i 

the number of terms the formula. By construction, the sets nw{k) and sw{l) 
are not empty and all terms are positive. In statistical physics and more specifi- 
cally in quasicrystal science, people are interested in thermodynamic quantities 
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(22) 



such as the configurational entropy (per tile): S = \n{Ta,b,c,d) / Nt where Nt 
is the number of tiles. With our polygonal boundary conditions, this quantity 
has a finite limit when goes to infinity provided the relative ratios of the 
side lengths also have a finite limit [1,2,10,12]. In the so-called "diagonal" case 
where all side lengths are equal, taking into account the number of tiles, the 
previous relation becomes S > S^^"" = (3/2) ln3 - 2 ln2 ~ 0.262 [12]. The ac- 
tual value of S is numerically known to be close to 0.36 [2]. The previous lower 
bound is manifestly loose and its improvement requires a better knowledge of 
the asymptotic behavior of the determinants in (12) at the large size limit. 

But the main advantage of our formula precisely lies on the fact that the previ- 
ous bound is weak: the formula realizes an exponential reduction of the number 
of terms as compared to a crude enumeration of tilings. Indeed, as it was just 
discussed in the previous paragraph, the number of terms grows exponentially 
hke exp(0.26A^T) wheres the number of tihngs grows like exp(0.36A^r)- Even if 
in practice we cannot compute numerically the number of tilings of octagons 
bigger than in table 1, the progress is already significant. Moreover there exists 
some hope to simplify our formula, at least partially, as in eq. (3). 

In addition, the formula brings a new insight into the structure of tiling sets: 
it emphasizes a natural decomposition of the sets into smaller disjoint subsets, 
the cardinality of which is simply given by evaluation of determinants. 
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